Abstract 

We present a search for the decays — + e+e^, M^A* , and 

e^/i^ in data collected at the T(4iS') resonance with the Babar 
detector at the SLAC B Factory. Using a data set of 111 fb~^, we find 
no evidence for a signal in any of the three channels investigated and 
set the following branching fraction upper limits at the 90% confidence 
level: B{B'^ e+e") < 6.1 x 10-^ B{B° M+zx") < 8.3 x 10-^ and 
B{B" e±/z=F) < 18 X 10-8. 
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We present a search for the decays e'^e~, — > /i'^ fi~ , and e^fi^ in data collected 

at the T(45') resonance with the BABAR detector at the SLAC B Factory. Using a data set of 
lllfb"^, we find no evidence for a signal in any of the three channels investigated and set the 
following branching fraction upper limits at the 90% confidence level: B{B° e+e") < 6.1 x 10"*, 
B{B° m^M") < 8.3 X 10"^ and B{B" e'^/i'^) < 18 x 10"^ 

PACS numbers: 13.20.He,14.40.Nd 



In the Standard Model (SM), rare B decays such as 
5° i^£^, where £ refers to e or /x, are expected 
to proceed through diagrams such as those shown in 
Fig. n (charge conjugate processes are included implicitly 
throughout). These decays are highly suppressed since 
they involve a 6 — > d transition and require an internal 
quark annihilation within the B meson. In addition, the 
decays are helicity suppressed by factors of {mi/mBY 
where is the mass of the lepton £ and tob is the mass 
of the B meson. i3° decays to leptons of two different fla- 
vors, violate lepton flavor conservation and are therefore 
forbidden in the SM, although permitted in extensions to 
the SM with non-zero neutrino mass . The SM expec- 
tations are given in Tabled 

Since these processes are highly suppressed in the 
SM, they are potentially sensitive probes of physics be- 
yond the SM. In the minimally supersymmetric Standard 
Model (MSSM) the branching fraction for these decays 
can be enhanced by orders of magnitude In par- 
ticular, for MSSM models with modified minimal flavor 
violation (MFV) and large tan/3 Q, the branching frac- 
tion can be increased by up to four orders of magnitude. 
Experimental bounds can restrict allowed regions of pa- 
rameter space, specifically the mass of the charged Higgs 
boson. In non-MSSM models with two Higgs doublets 
and natural flavor conservation at large tan /3, an increase 
in the branching fraction of several orders of magnitude 
is expected 01 ■ B° ^ decays are also allowed in 

specific models containing leptoquarks [5||and supersym- 
metric (SUSY) models without R-parity [6j. The branch- 
ing fractions for the flavor violating channels ^t^J 
(i j) are expected to be exceedingly small but can be 
enhanced by leptoquarks or _R-parity violating operators 
in SUSY models. 

To date, B° V^l^ decays have not been observed. 
As shown in Tabled experimental limits are approaching 
a level of sensitivity that will restrict the allowed param- 
eter space of models that produce i?° — > branching 
fraction enhancements of a few orders of magnitude with 
respect to the SM rates. 

The data used in these analyses were collected with 
the BABAR detector at the PEP-II e+e~ storage ring 
and correspond to an integrated luminosity of 111 fb~^ 
accumulated at the T(\S) resonance ( "on-resonance" ) 
and 11.9 fb~^ accumulated at a center-of-mass (cm.) 
energy about 40MeV below the T(\S) resonance ("off- 




FIG. 1: Representative Feynman diagrams for B'^ £^£ in 
the Standard Model. 



TABLE I: The expected branching fractions in the Standard 
Model ^ and the current best upper limits (U.L.) at the 90% 
C.L. 



Decay 
Mode 


SM 
Expectation 


CLEO \8] 
g.lfb"^ 


Belle [2] 
78 fb"^ 


CDF [12] 
0.17fb"i 


e+e- 


1.9 X 10"^^ 


8.3 X lO"'^ 


1.9 X lO"'^ 






8.0 X 10"" 


6.1 X lO"'^ 


1.6 X 10"'^ 


1.5 X 10"'^ 


e^fi^ 




15 X lO"'^ 


1.7 X lO"'^ 





resonance"). The latter sample is used for non- resonant 
QQ {q = fJ-^ s, and c) background studies. The collider 
is operated with asymmetric beam energies, producing a 
boost (/37 — 0.55) of the T(4iS') along the collision axis. 

The BABAR detector is optimized for the asymmetric 
beam configuration at PEP-II and is described in de- 
tail in 11]. The 1.5-T superconducting solenoidal mag- 
net, whose cylindrical volume is 1.4 m in radius and 
3 m long, contains a charged-particle tracking system, 
a Cherenkov detector dedicated to charged-particle iden- 
tification, and central and forward electromagnetic Csl 
calorimeters (EMC). The segmented flux return, includ- 
ing endcaps, is instrumented with resistive plate cham- 
bers for muon and K'^^ identification. 

The presence of two charged high-momentum leptons 
provides a very clean signature for the three decay modes 
under consideration. We require two oppositely-charged 
high- momentum leptons {i.e. \p*^\ ^ ttibI'^, where p| is 
the cm. momentum of lepton ^) from a common ver- 
tex consistent with the decay of a B'^ meson. Since the 
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signal events contain two mesons and no additional 
particles, the total energy of each B° in the cm. must 
be equal to half of the total beam energy. We define 



rriES 



AE = Y^^mj + ivlY-E, 



(1) 

beam ' (2) 



where -Ebcam (^^ ) beam energy in the cm. 

frame, p* is the momentum of lepton i in the cm. frame, 
and rrii is the mass of lepton i. In Eq. -E^cam used 



as opposed to E'^ because i?boam known with much 
greater precision. For correctly reconstructed mesons. 
Toes peaks at the mass of the meson with a resolution 
of about 2.5 MeV/c^ and AE peaks near zero. 

To reduce background from lepton misidentification, 
we require the leptons to satisfy stringent electron and 
muon identification criteria |l2l |. The electron iden- 
tification efficiency is greater than 93% with a mis- 
identification rate of less than 0.3%. The muon identifi- 
cation efficiency ranges from (55 — 70)% (depending on 
run period) with a mis- identification rate of 3%. Electron 
energy lost through bremsstrahlung is partially recovered 
by adding the energy of photons that lie within a 3 degree 
cone about the electron direction. 

Suppression of background from non-resonant qq pro- 
duction is provided by a series of topological require- 
ments. In particular, we require |cos0t| < 0.8, where 
6t is the angle in the cm. frame between the thrust axis 
of the particles that form the reconstructed B^ candidate 
and the thrust axis of the remaining tracks and neutral 
clusters in the event. In addition we employ cuts on the 
invariant mass of the "Rest Of the Event" (all tracks not 
associated with the B^ candidate where all non-leptonic 
tracks are assumed to be pions) of torqe > 0.5 GeV 
and on the second normalized Fox- Wolfram moment of 
i?2 < 0.8 01 . We also cut on the total multiplicity 
of both charged tracks and neutral particles by means 
of the variable iVmuit defined as A^muit = -^trk + 
where iVtrk is the total number of tracks in the event 
and is the number of photons found with an energy 
E^ > 80 MeV. We require TV^uit > 5.5 for the ee and 
e/i channels and N^^it > 5.0 for the fifi channel. This 
variable is especially useful in the rejection of radiative 
Bhabha events. We also require that the total energy in 
the EMC (E'emc) be less than 11 GeV. This cut is ef- 
fective in reducing background from QED e"^e~ events, 
including radiative Bhabhas with many conversions. 

Four of the selection criteria given above {\ cos 9t\, 
™ROE, -^muit, and -Eemc) were simultaneously optimized 
for the best upper limit on B{B'~' — > i^£~) where the as- 
sumed number of observed events is determined from a 
Poisson distribution with the mean equal to the expected 
background. Sideband data are compared with signal 



mO.3 




FIG. 2: Distributions of signal MC (hatched) and sideband 
data (points) for the e^e~ channel after the initial selection 
cuts for (a) | cos0t|, (b) rriROE, (c) R2, and (d) A^'muit- Arrows 
indicate final cut values. All distributions are normalized to 
unity. 



Monte Carlo (MC) for the e+e^ channel for four of these 
variables, |cos6't|, torqe, R2, and A'muit in Fig. [5] 

The B'^ candidates are selected by simultane- 

ous requirements on the energy difference AE and the 
energy-substituted mass toes- For the B^ de- 
cay mode, the size of this "signal box" is chosen to be 
[+2, —2]a of the toes distribution and [-1-2, — 2](t for the 
AE. In the cases of the B° e+e" and B° e^^i^ de- 
cay modes, the signal box sizes in toes are also [-1-2, — 2](t 
but in AE are relaxed to [4-2, — 3]cr and [-1-2, — 2.5]o', re- 
spectively, to accommodate the tail in the distribution re- 
sulting from uncorrected bremsstrahlung and final state 
radiation. The resolution in toes is obtained from a fit 
to a Gaussian distribution, whereas the resolution in AE 
is obtained from a fit to an empirical function [14] that 
gives a good description of this tail. 

We estimate the background level in the signal box 
from the data sidebands assuming that it is described 



TABLE 11: Definition of the three different sideband boxes 
used for all three decay modes. 



Sideband Box 


span in AE [GeV] 


span in rriES [GeV/c'^] 


Grand Sideband 
Upper rriES 
Lower jtibs 


(-0.40,0.40) 
(0.20,0.50) 
(-0.50, -0.20) 


(5.20,5.26) 
(5.20,5.29) 
(5.20,5.29) 
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TABLE III: Summary of the analyses where A^obs and TV^p 
are the observed and expected number of events in the signal 
box, e is the efficiency, and Bul(B° — > £'^£~) is the upper 
limit on the branching fraction at the 90% C.L. Systematic 
uncertainties on N^^p and e are given. 



Decay Mode A'obs 


' cxp 


s[%] 


Bui^iB" ^ £+£-) 


e+e" 


0.71 ±0.31 


21.8 ± 1.2 


6.1 X IQ-* 





0.72 ± 0.26 


15.9 ± 1.1 


8.3 X IQ-* 


e±/i=F 2 


1.29 ± 0.44 


18.1 ± 1.2 


18 X 10"* 



by the ARGUS function \l$ in mES and an exponen- 
tial function in AE. We use these parameterizations to 
extrapolate the background level found in the sidebands 
into the signal box. As indicated in Table ITll three differ- 
ent sideband boxes are used. The grand sideband box is 
used to estimate the functional form of the AE distribu- 
tion. The upper and lower tues sideband boxes are used 
to estimate the functional form of the tties distribution. 
Peaking backgrounds from misidentified two-body B de- 
cay modes were estimated using an MC sample equiva- 
lent to more than 20 times the data luminosity and found 
to be negligible. The total background expectations and 
signal efficiencies are given in Table Hill 

The systematic uncertainties on the signal efficiency e, 
the number of B'^ mesons produced in the data, and the 
background estimate are incorporated into the determi- 
nation of the upper limit on B{B'^ — ^ i^i^). Since the 
signal efficiency is determined from MC simulation only, 
differences between data and the simulation would result 
in an error in our normafization. To estimate this un- 
certainty we perform comparisons of data and MC using 
high statistics control samples that have similar charac- 
teristics to our B'^ t^l~ signal. The optimal con- 
trol samples are B^ — > J/ipK^^ with J/%Ij ^ e+e~ for 
e+e" and J/il) ^ for B^ M^M > respec- 

tively. Since there exists no appropriate control sample 
for the e^fi^ mode, we use the larger of the system- 
atic errors derived from either the ee or /i/x modes. In 
performing these comparisons we found a substantial un- 
certainty on the signal efficiency to be due to differences 
between data and the MC simulation in the mean and 
resolutions of various quantities, depending on the chan- 
nel. For the electron channels the dominant quantities 
are AE and mROE whereas for the muon channels they 
are | cos 9t \ , N^^n , and m^oE ■ When combined with 
the uncertainties on tracking efficiency of 2.6% and that 
for particle identification (1.0% per electron, 3.0% per 
muon) , the total systematic uncertainty on the efficiency 
is estimated to be 5.7%, 7.1%, and 6.8% for the ee, fifi, 
and e/i modes respectively. 

The background estimate is obtained from a fit to side- 
band data, so the primary uncertainty is due to fluctua- 
tions in the fit procedure as events fall in or out of the 



sideband box. We have studied the stability of the fit and 
the background estimate when adding or removing events 
from the thes and AE histograms. We find that the fit is 
unbiased and stable to a level significantly less than the 
statistical uncertainty on the background estimate. 

As shown in Fig. O and Table IIIII when the contents 
of the signal box were revealed, 0, 0, and 2 events were 
found in the ee, /i/i, and e^ channels respectively. As 
can be seen in Table IIIII the numbers of events found in 
the signal boxes are compatible with the expected back- 
ground for each mode. 

An upper limit on the branching fraction is computed 
using 

B{B°^t+£-)= , (3) 

(Abo + Ngo ) ■ e 

where A[/L(Aobs) is the Poisson 90% U.L. on the num- 
ber of events assuming Aobs events have been observed, 
NBo{Ngo) is the number of B^{B^) mesons produced 
in the data, and e is the signal efficiency. We have 
Nbo -\- Ngo = uudcr the assumption of equal pro- 

duction of B^B" and B+B^ in T(4S') decays. For our 
data set N^g = (122.5 ± 1.2) x 10^. 

We follow the technique of jl^] in order to account for 
the presence of background and to include our systematic 
uncertainties in the determination of the upper limit. As 
summarized in Table IIIII the resulting upper limits at 
the 90% C.L. for B{B° e+e"), B{B° H+H'), and 
B{B° e±AiT) are 6.1 x 10-^ 8.3 x 10-^ and 18 x 10"^ 
respectively. 

These bounds are stringent enough to place interest- 
ing constraints on popular models. For example, for the 
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FIG. 3: Distribution of events in mBS and AE for B" — > e+e 
(top left), B° ^i+fj,- (top right), and e+^t+ (bottom). 
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MSSM (MFV) models, the relation between B{B° 
fi^fi~) and the mass of the charged Higgs boson mn 
is given as a function of tan/3 in j^. We find that for 
tan/? = 60, ruH > 138 GeV (90% C.L.). 
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